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Abstract 

Weyl theory for Dirac systems with rectangular matrix potentials 
is non-classical. The corresponding Weyl functions are rectangular 
matrix functions. Furthermore, they are non-expansive in the upper 
semi-plane. Inverse problems are treated for such Weyl functions, and 
some results are new even for the square Weyl functions. High energy 
asymptotics of Weyl functions and Borg-Marchenko type uniqueness 
results are derived too. 

MSC(2010): 34B20, 34L40. 

Keywords: Weyl function, Weyl theory, Dirac system, rectangular matrix 
potential, inverse problem. 

1 Introduction 

The self-adjoint Dirac-type (also called Dirac, ZS or AKNS) system 

-^y{x,z) =i{zj + jV{x))y{x,z) (x > 0), (1.1) 
ax 



where 



J 



Im, 

-^rn.o 



V 



V 
V* 



:i.2) 



/mj. is the rrifc x nik identity matrix and v{x) is an mi x m2 matrix function, is a 
classical matrix differential equation. We always assume that v is measurable 
and, moreover, locally summable, that is, summable on all the finite intervals. 
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Dirac-type systems are very well-known in mathematics and applications, 
especially in mathematical physics (see, e.g., books [7,8,26,28,43], recent 
publications [4-6,9,10,16,17,44], and numerous references therein). System 
(1.1) is used, in particular, in the study of transmission lines and acoustic 
problems [48]. The most interesting applications are, however, caused by 
the fact that system (1.1) is an auxiliary linear system for many important 
nonlinear integrable wave equations and as such it was studied, for instance, 
in [1-3,13,18,23,24,36,49]. (The name ZS-AKNS system is connected with 
these applications.) Nonlinear Schrodinger equations, modified Korteweg de 
Vries equations, and second harmonic generation model, which describe vari- 
ous wave processes (including, e.g., nonlinear wave processes in water and in 
waveguides, in nonlinear optics, photoconductivity, and on silicon surfaces), 
are only some of the well-known examples. The evolution of Weyl func- 
tion for these systems is described in terms of Mobius transformations (see, 
e.g., [32,36,37,41-43] and references therein), which is one of the fruitful ap- 
proaches to study interesting initial-boundary value problems for integrable 
nonlinear equations. 

The Weyl and spectral theory of self- adjoint Dirac systems, where mi = 
777.2, was treated, for instance, in [6,9,25,28,35,43] (see also various references 
therein). The "non-classical" Weyl theory for the equally important case 
nil 7^ ^2, which appears in the study of coupled, multicomponent, and 
matrix nonlinear equations, is the subject of this paper. 

The mi x m2 matrix function v{x) from (1.2) is called the potential of 
system (1.1). We put mi + m2 =: m. The fundamental solution of system 
(1.1) is denoted by u{x, z), and this solution is normalized by the condition 

U{0,z)^lm. (1.3) 

The next Section 2 is dedicated to representation of the fundamental 
solution. In Section 3 we follow [14] to introduce the Weyl function, and 
study the high energy asymptotics of this Weyl function. The solution of the 
inverse problem to recover the potential (and system) from the Weyl function 
is given in Section 4. Borg-Marchenko-type uniqueness results are contained 
in that section as well. Finally, Section 5 is dedicated to conditions for an 
analytic matrix function to be the Weyl function of some Dirac-type system. 
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The results from this section were not pubhshed before even for the case that 
mi = m2, though some of them appeared for that case in the thesis [34]. 

As usual, C stands for the complex plain, C+ for the open upper semi- 
plane, and Im stands for image. If a G C, then a is its complex conjugate. 
The notation L^^xmilO) ^) '^i^ be used to denote the space of 1712 x mi 
matrix-functions with entries belonging to L^(0, 00). We use H to denote the 
closure of the space H, I to denote the identity operator, and -B(Hi, H2) to 
denote the class of bounded operators acting from Hi to H2. We write -B(Hi) 
if Hi = H2, and -B^[0, I] will be used to denote a class of functions (or matrix 
functions), whose derivatives are bounded on [0, I]. An m2 x mi matrix a is 
said to be non-expansive, if a*a < (or, equivalently, if aa* < 1^.2)- 

2 Representation of fundamental solution 

The results of this section can be formulated for Dirac system on a fixed final 
interval [0, /]. We assume that v is bounded on this interval and put 

P{x)^[lm, 0]U{X,0), ^{X)^[0 Im2]u{x,0). (2.1) 

It follows from sup^.^; II^H < oo and from (1.1) that 

sup||7'(a;)|| < oo, 7' := -^7. (2.2) 
x<i 

Moreover, by (1.3) and (3.2) we get u{x,0)*ju{x,0) = j = u{x,0)ju{x,Oy . 
Therefore, (2.1) implies 

/3j/3* = V, 7jY = -/^^, /3n* = 0. (2.3) 

Next, we need the following similarity result for the Volterra operator 

r F{x)G{t) -dt, Ke B{Ll^{0, I)), (2.4) 
^0 

where F{x) is an mi x m matrix function, G{t) is an m x mi matrix function, 
and 

iF{x)G{x) = I^,. (2.5) 
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Proposition 2.1 [38] Let F and G he houndedly differentiable and let (2.5) 
hold. Then we have 

K = EAE-\ A:=-ir-dt, A, E, e B{LI^{0, I)) , (2.6) 

^0 

where K is given by (2.4) and E is a triangular operator of the form 

{Ef){x) = p{x)f{x) + J E{x, t)f{t)dt, = iF'Gp, det p(0) ^ 0. 

(2.7) 

Moreover, the operators E^^ map functions with hounded derivatives into 
functions with bounded derivatives. 

The proposition above is a particular case of Theorem 1 [38] (see also a later 
paper [5] on the case of continuously differentiable F and G) . 
Set 

F = 7, G-ijT*, 7J7* = -/™„ 7e5^[0, /], (2.8) 

where -B^ stands for the class of houndedly differentiable matrix functions. 
Clearly F and G in (2.8) satisfy conditions of Proposition 2.1. Separate 
7(a;) into two blocks 7 = [71 72], where 71, 72 are m2 x mi and m2 x m2, 
respectively, matrix functions. Without loss of generality one can choose E 
so that 

E-^^2 = Im„ (2.9) 
where E"^ is applied to 72 columnwise. 

Proposition 2.2 Let K be given by (2.4), where F and G satisfy (2.8), 
and let E be a similarity operator from Proposition 2.1. Introduce Eq e 
B(^L'^^{0, /)) by the equalities 

{Eof){x) = p(0)-S2(0)/(x) + rEoix - t)f{t)dt, Eo{x) := {E-'^^Yix). 

Jo 

(2.10) 

Then, the operator E := EEq is another similarity operator from Proposition 
2.1, which satisfies the additional condition (2.9). 
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Proof. The proof of the proposition is similar to the case mi = m2 (see, 
e.g., [38, pp. 103, 104]). Indeed, the next identity can be easily shown directly 
(and follows also from the fact that Eq is a lower triangular convolution 
operator) : 



^EqA, (A^-i f -dt). (2.11) 
Jo 



AEo 

'0 



Furthermore, because of the third relation in (2.8) we have det72(0) ^ 0, and 
so £^0 is invertible. Hence, equalities (2.11) and K = EAE~^ imply (2.6): 

K = EAE-^, E := EEq. (2.12) 

Formula (2.10) leads us also to the equality 

(Eolm,) (x) = p(0)-S2(0) + / Eo{x - t)dt = p(0)-S2(0) + / E^{t)dt 

Jo Jo 

= p(0)-S2(0) + {E-'^2){x) - {E-'^2){0). (2.13) 
Recalling that E~^'j2 G -S^[0, /] we obtain a representation: 

{E-'j2)ix) = p(0)-S2(0) +i(^(^-S2)')(^). (2.14) 
Using (2.14) we rewrite (2.13) as 

{EoIm,){x) = {E-'^2){x), (2.15) 

and (2.9) follows. It remains to show that E^^ maps -B^[0, /] into B^[0, /]. 
First note that the integral operators Eq, Eq^ have bounded kernels and 
map bounded functions into bounded. In particular, for Eq^ it follows from 

a series representation of the operator of the form ^/ + k{x,t) ■ dt^ 
(where the kernel k is bounded). Now, taking into account that (similar to 
(2.14)) any / G B^[0, I] admits representation / = /(O) + iAf, we see that 
formulas (2.11) and (2.15) yield the fact that Eq maps ^^[O, /] into B^[0, I]. 
The identity AE^^ — Eq^A is immediate from (2.11). Thus, to prove that 
Eq^ maps Si[0, l\ into ^^[0, I] we need only to show that Eq ^I^^ e -B^[0, I], 
which can be derived from (2.14) and (2.15): 

E-'I^, = (/^, -iEo"'^(^"S)')72(0)'V(0) e B% I]. 

□ 
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Reiricirk 2.3 The kernels of the operators -E , which are constructed in 
[38], as well as the kernels of the operators E^^ from the proof of Proposition 
2.2 are bounded. Therefore, without loss of generality we always assume 
further that the kernels of E^^ are bounded. 

The next lemma easily follows from Proposition 2.2 and will be used to 
construct fundamental solution. 

Lemma 2.4 Let 7 be an m2 x m matrix function, which satisfies the last 
two relations in (2.8), and set 

S:=E-'{E*)-\ n:=[$i$2], G i?(C™^ L^^(0, /)); (2.16) 

($i/)(a;) = $i(x)/, $i(x)-(£;-S)(^); ^2f^ImJ = f; (2.17) 

where E is constructed (for the given ^) in Proposition 2.2. Then A, S, and 
n form an S-node, that is (see [39,40,43]), the operator identity 

AS - SA* = ilijW (2.18) 

holds. Furthermore, we have 

00 

Y,^m{{A*YS-^Tl)^Ll^{Q,l). (2.19) 

Proof. Because of (2.4), (2.6) and (2.8), we get 

EAE-^ - {E-yA*E* = K-K* = i-^{x)j [ ^{t)* ■ dt. (2.20) 

Jo 

Formulas (2.9), (2.16) and (2.17) lead us to the equality 

n/=(£;-S)(x)/. (2.21) 

Now, the operator identity (2.18) follows from (2.20), (2.21), and the first 
equality in (2.16). 

To prove (2.19) we will show that 

N N N 

J]lm (^{A*yS-^U^ D ^Im {S-^A'U) = 5'^ ^ Im {A'U). (2.22) 

i=0 i=0 i=0 
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For that purpose we rewrite (2.18) as S'^A = A* S'^ +iS-^UjU* S'^ . Hence, 
for Ni, N2 > we obtain 

Im (^{A*f"^'s-'A''m^ (2.23) 
N1+N2 

+ J] Im (^{A*yS-'uj Dim (^{A*f'S-'A'''+'uy 

1=0 

Using (2.23), we derive (2.22) by induction. In view of (2.22), it suffices to 
show that 

00 

J]lm(A^n)=L^^(O,0, (2.24) 

i=0 

which, in its turn, follows from (2.17). □ 

Remark 2.5 Given an S-node (2.18), we introduce a transfer matrix func- 
tion in Lev Sakhnovich form {see [39,40,43]) '■ 

WA{r, z) := + izjU*S;\l - ^^)"^P,.n, < r < Z, (2.25) 

where I is the identity operator; A^., Sr G B(^L'^^{0, r)), 

Ar := PrAP;, Sr := PrSP;-, (2.26) 

A is given by (2.6), the operators S and U are given by (2.16) and (2.17), 
and the operator Pr is an orthoprojector from L'^^{0, I) on L'^^{0, r) such 
that 

{Pj){x) = f{x) {0<x<r), feLl^iO,l). (2.27) 

Since P^A = P.f.AP*P.f., it follows from (2.18) that the operators A^, S^, and 
Pril form an S-node too, that is, the operator identities 

ArSr - SrA; = iP^UjU* P^ (2.28) 

hold. 
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Now, in a way similar to [31,32] the fundamental solution w of the system 

z) = izj-f{xy-f{x)w{x, z), w(0, z) = (2.29) 

is constructed. 

Theorem 2.6 Let 7 be an m2 x m matrix function, which satisfies the last 
two relations in (2.8). Then, the fundamental solution w given by (2.29) 
admits representation 

w{r,z) = WA{r,z), (2.30) 
where WA{r, z) is defined in Remark 2.5. 

Proof. The statement of the theorem follows from Continual factoriza- 
tion theorem (see [43, p. 40]). More precisely, our statement follows from 
a corollary of the Continual factorization theorem, namely, from Theorem 
1.2 [43, p. 42]. Using Lemma 2.4 we easily check that the conditions of 
Theorem 1.2 [43, p. 42] are fulfilled. Therefore, if n*5'~^P^n is boundedly 
differentiable, we have 

-^WA{r, z) = \zjH{r)wA{r, z), lim ^^(r, z) = (2.31) 
ar r-^+o 

H(r) := ^{U*S;:'PrU), (2.32) 

where wa is given by (2.25). Since E"^^ are lower triangular operators, we 
see that 

PrEP^Pr = PrE, {E-^) * P^ = P^ Pr{E-^) * P,* . (2.33) 

Hence, formulas (2.16) and (2.26) lead us to 

S;^ = E;Er, Er := PrEP;. (2.34) 
Therefore, taking into account (2.21), we rewrite (2.32) as 

H{r) = 7(r)*7(r). (2.35) 
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Formulas (2.29), (2.31) and (2.35) imply (2.30). □ 

Now, consider again the case of Dirac system. Because of (2.1) and (2.3) we 
obtain 

u{x,0)jj{x)*^{x)u{x,0)'^ — — 

Hence, direct calculation shows that the following corollary of Theorem 2.6 
is true. 

Corollary 2.7 Let u{x, z) be the fundamental solution of a Dirac system 
with a bounded potential v and let 7 be given by (2.1). Then u{x,z) admits 
representation 

u{x, z) = e^'^'uix, 0)w{x, 2z), (2.37) 

where w has the form (2.30) and the S-node generating the transfer matrix 
function wa is recovered from 7 in Lemma 2.4- 

Remark 2.8 For the case that j is given by (2.1), it follows from (1.1) and 
(2.3) that 

y{x)n{x)* = -i [v{x)* 0] u{x, 0)n{x)* = -w{x)*/3{x)n{x)* = 0. (2.38) 

Thus, from (2.8) and (2.38) we see that p in (2.7) is a constant matrix. 
Therefore, since 72(0) = /m27 equality (2.21) implies that p{x) = Im2, o-nd 
formula (2.7) can be rewritten in the form 

{Ef){x)^f{x)+ r E{x,t)f{t)dt, EeB{Ll^{0,l)). (2.39) 
^0 

RecaUing that det72(x) 7^ we can rewrite (2.38) as 

7; = 7;(72"S)*- (2.40) 
Using (2.40), we recover 72 and 71 = 72(72"^7i) from 7^^7i. Indeed, we have 

(72"Si)' = -72"S272"Si + 72"S1> i-e-, 7i = 72(72"Si) + 72(72"Si)'- 





, 



(2.36) 
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Because of (2.1), (2.40) and formula above we get 

72(/™, - (72"Sl)(72~Sl)*) = 72(72"Sl)'(72"Sl)*, 72(0) = J™,. 

Therefore, taking into account that 1^2 ~ (72'^7i)(72'"^7i)* — 72'"^(72'"^)* 
invertible, we obtain 

72 = 72(72"Si)'(72"Si)*(/m. - (72"Sl)(72"Sl)*)~\ 72(0) = (2.41) 

Remcirk 2.9 If the conditions of Remark 2.8 hold, then 72, 7 and Hamil- 

tonian H = 7*7 are consecutively recovered from 7^^7i via the differential 
equation (2.41). Furthermore, from the third relation in (2.8) we see that 
(72^^71) (72^^71)* < Im2- Thus, the matrix function 7^^7i is a continuous 
analog of Schur coefficients for the canonical system (2.29). (Compare with 
Remark 3.1 from [17] on the canonical systems, which were treated there.) 

3 Weyl function: high energy asymptotics 

To define Weyl functions, we introduce a class of nonsingular m x mi ma- 
trix functions V{z) with property- j, which are an immediate analog of the 
classical pairs of parameter matrix functions. Namely, the matrix functions 
Viz) are meromorphic in C+ and satisfy (excluding, possibly, a discrete set 
of points) the following relations 

V{zyV{z) > 0, V{z)*jV{z) > {ze C+). (3.1) 

It is immediate from (1.1) that 

— {u{x, z)*ju{x, z)) — \{z — z)u{x, z)*u{x, 2;) < 0, z E C+. (3.2) 

Relations (3.1) and (3.2) imply 

det([V 0]u{x,z)-'V{z)') 7^0, (3.3) 

Definition 3.1 The setAf{x,z) of Mohius transformations is the set of val- 
ues at X, z of matrix functions 

ipix,Z,V)=[0 Im2]uix,z)~'V{z)(^[lrr.^ O] uix, z)-'V (z)^ \ (3.4) 

where V{z) are nonsingular matrix functions with property- j. 
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We can rewrite (3.4) in an equivalent form, which will be used later on 

u{x, z)-^V{z) ( [Irr,, O] u{x, z)-'V{z)) '\ (3.5) 



I mi 



if{x,z,V)^ 

Proposition 3.2 [14] Let Dirac system (1.1) on [0, oo) be given and assume 
that V is locally summahle. Then the sets J^{x, z) are well-defined. There is 
a unique matrix function ip{z) in C+ such that 

^{z) = fl M{x,z). (3.6) 

x<oo 

This function is analytic and non-expansive. 

In view of Proposition 3.2 we define the Weyl function of Dirac system similar 
to the canonical system case [43]. 

Definition 3.3 [14] The Weyl-Titchmarsh (or simply Weyl) function of 
Dirac system (1.1) on [0, oo), where potential v is locally summahle, is the 
function </? given by (3.6). 

By Proposition 3.2 the Weyl-Titchmarsh function always exists. Clearly, it 
is unique. 

Corollary 3.4 [I4] Let the conditions of Proposition 3.2 hold. Then the 
Weyl function is the unique function, which satisfies the inequality 



POO 

/ [Imi ip{z)*]u{x,zyu{x,z) 

Jo 



dx < 00. (3.7) 



Remark 3.5 In view of Corollary 3.4, inequality (3.7) can be used as an 
equivalent definition of the Weyl function. Definition of the form (3.7) is a 
more classical one and deals with solutions of (1.1) which belong to L^(0, 00). 
Compare with definitions of Weyl-Titchmarsh or M -functions for discrete 
and continuous systems in [10, 28, 29, 32, 33, 43, 46, 47] (see also references 
therein). 



11 



Important works by F. Gesztesy and B. Simon [19,20,45] gave rise to a 
whole series of interesting papers and results on the high energy asymptotics 
of Weyl functions and local Borg-Marchenko-type uniqueness theorems (see, 
e.g., [9,11,12,17,27,35] and references therein). Here we generalize the high 
energy asymptotics result from [35] for the case that Dirac system (1.1) has a 
rectangular mi x 1712 potential v, where mi is not necessarily equal to 1712. For 
that we recall first that 5" > and $1 is boundedly differentiable in Lemma 
2.4. Therefore, using Theorem 2.5 from [15] we get the statement below. 

Theorem 3.6 Let U = [$1 $2] be constructed in Lemma 2.4- Then there 
is a unique solution S G B(^L'^^{0,1)^ of the operator identity (2.18), this S 
is strictly positive {i.e., S > 0) and is defined by the equalities 



Continuous operator kernels of the form above (with a possible jump at 
X — t) were considered in Section 2.4 from [21], where they were called "close 
to displacement kernels" (see also references therein). 

Now, we will apply the (S-node scheme, which was used in [32] for the 
skew-self- adjoint case, to derive the high energy asymptotics of the Weyl 
function of Dirac system with a rectangular potential. 

Theorem 3.7 Assume that ip e J\f{l,z), where J\f{l,z) is defined in Def- 
inition 3.1 and the potential v of the corresponding Dirac system (1.1) is 
hounded on [0, I]. Then {uniformly with respect to '^{z)) we have 



{Sf){x) = (/^, - <l>i(0)$i(0)*)/(a;) - / s{x,t)f{t)dt, 



(3.8) 




X > t; 



t> X. 



(p{z) — 2iz / e 








— >■ 00. 



(3.9) 



Proof. To prove the theorem, we consider the matrix function 



U{z) - [V V{zr] {j - wa{1, 2zrjWA{l, 2z)) 



(3.10) 
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Because of (2.30) and (2.37) we have 

U{z) = V - ^{zr^{z) - e"(^-^) [V ^{zY] u{l, zYMl, , 



-'mi 



Taking into account (3.5), we rewrite the formula above as 
U{z) = Im, - ^{zT^iz) - e«(^-^) (3.11) 
x(([7^, 0]u{x,z)-'V{z))-yV{zrjV{z){[l^^ 0]u{x,z)-'V{z))-\ 

From (3.1) and (3.11) we see that 

U{z)<Im,. (3.12) 
It easily follows from (2.18) and (2.25) (see, e.g., [17,40]) that 

wa{1, zfiwAil, z)=j + i{z - z)I{*{I - zA*Y^S^\l - zAY^n. (3.13) 
Now, formulas (3.10), (3.12) and (3.13) imply 



2i(^ - z) [Im, ^{z)*] n*(/ - 2zA*)-^S;\l - 2zA)-^Ii 



Since S is positive and bounded, inequality (3.14) yields 



-'mi 



(/ - 2zA)-^Il 



ip{z) 



< C/V^ for some C > 0. 



mi • 

(3.14) 
(3.15) 



We easily check directly (see also these formulas in the works on the case 
mi = 777.2) that 



mi -2zA)-^f ^ I e^'^''-^^'f{x)dx, 
Jo 

mi - 2zA)-'^2 = ^(e"'"' - 1)^- 



(3.16) 
(3.17) 



Because of (3.15)-(3.17) (and after applying — i$2 to the operator in the 
left-hand side part of (3.15)), we get 



2z 



(e-2"^-l)^(z) =ie- 



-2ilz 



I e^^^^^i 
Jo 



{x)dx + O 



(3.18) 
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Since (/? is non-expansive, we see from (3.18) that (3.9) holds. □ 

Now, consider a potential which is locally bounded, that is, bounded on 
all the finite intervals [0, /]. The following integral representation is essential 
in interpolation and inverse problems. 

Corollary 3.8 Let ip he the Weyl fuction of Dirac system (1.1) on [0, oo), 
where the potential v is locally bounded. Then we have 

/•oo 

ip{z) = 2iz / e^^^^$i(x)dx, ^{z) > 0. (3.19) 
Jo 

Proof. Since (p is analytic and non-expansive in C_|_, for any £ > it admits 
(see, e.g., [30, Theorem V]) a representation 

/■oo 

(fi{z) = 2iz / e^^^^$(x)dx, ^{z) > £ > 0, (3.20) 
Jo 

where e"^^^$(a;) e L'^^^^^{0, oo). Because of (3.9) and (3.20) we get 

Q{z) : = [ e2'(^-')^($i(a;) - ^{x))dx 
Jo 

/oo 
e^^i^-l>^x)dx + 0{l/^/^z)). (3.21) 

From (3.21) we see that Q{z) is bounded in some semi-plane '^{z) >r)o>0. 
Clearly, Q{z) is bounded also in the semi-plane Q{z) < rjo- Since Q is analytic 
and bounded in C and tends to zero on some rays, we have 

Q{z)= I e2'(^-')"($i(a;)-$(a;))(ia; = 0. (3.22) 
Jo 

It follows from (3.22) that $i(x) = $(x) on all the finite intervals [0, l\. 

Hence, (3.20) implies (3.19). □ 

Reiricirk 3.9 Since $i = we get that $i(a;) does not depend on I fori > x. 
Compare this with the proof of Proposition 4-1 in [17], where the fact that 
E{x, t) (and so ^i) does not depend on I follows from the uniqueness of the 
factorizations of operators S^^. See also Section 3 in [5] on the uniqueness of 
the accelerant. Furthermore, since $i = $ the proof of Corollary 3.8 implies 
also that e~^^^i{x) e L'^^xmii^^ f^''" ^''^V £ > 0. 
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Remcirk 3.10 From (3.19) we see that ^'^ is a Dirac system analog of A- 
amplitude, which was studied in [20, 45]. On the other hand ^[ is closely 
related to the so called accelerant, which appears for the case that mi = m2 
in papers by M. Krein (see, e.g., [5,25,26]). See also Remark 2.3 in [17], 
where $^ is discussed for the case that mi = m2. 

4 Inverse problem and B org- Mar chenko- type 
uniqueness theorem 

Taking into account Plancherel Theorem and Remark 3.9, we apply inverse 
Fourier transform to formula (3.19) and get 

t.(|)^W.._£e--||±M,5, „>o. (4.1) 

Here l.i.m. stands for the entrywise limit in the norm of iv^(0, r), < r < oo. 
(Note that if we put additionally $i(x) = for x < 0, equality (4.1) holds 
for l.i.m. as the entrywise limit in L^(— r, r).) Thus, operators S and 11 are 
recovered from (p. 

Since Hamiltonian H = 7*7 is recovered from S and 11 via formula (2.32), 
we recover also 7. Indeed, first we recover the Schur coefficient (see Remark 
2.9 for the motivation of the term "Schur coefficient"): 



[0 Im,]H 







[0 Im,]H 



(72*72) '72*71 = 72 S- (4.2) 



Next, we recover 72 from 7^ ^71 using (2.41), and finally, we easily recover 7 
from 72 and 7^ ^71. 

To recover /3 from 7, we separate /3 into two blocks /3 — [/3i ^2] , where 
f3k {k — 1,2) is an mi x m^ matrix function. We put 

^=[V 7i*(72*)^']- (4.3) 
Because of (2.3) and (4.3), we have /3jY — ^31* — 0, and so 

P{x) = Pi{x)^{x). (4.4) 
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It follows from (1.1) and (2.1) that 



(4.5) 



which implies 



P'jP* = 0, 



(4.6) 



Formula (4.4) and the first relations in (2.3) and (4.6) lead us to 




(4.7) 



-1 



(4.8) 



Relations (4.1)-(4.6) and (4.8) give us a procedure to construct the solution 
of the inverse problem. 

Theorem 4.1 Let he the Weyl function of Dirac system (1.1) on [0, oo), 
where the potential v is locally bounded. Then v can be uniquely recovered 
from (fi via the formula 



Here /3 is recovered from 7 using (4.3), (4.4) and (4.8); 7 is recovered from 
the Hamiltonian H using (4.2) and equation (2.41); the Hamiltonian is given 
by (2.32) andU and S in (2.32) are expressed via ^i{x) in formulas (2.17) 
and (3.8). Finally, ^i{x) is recovered from (p using (4.1). 

There is another way to recover ^ and 7: 

Remark 4.2 We can recover (5 directly from H and S following the propo- 
sition below, and then recover 7 from (3 in the same way that (3 is recovered 
from 7. 

Proposition 4.3 Let Dirac system (1.1) on [0, 00) he given. Assume that 
V is locally bounded and H and S are operators constructed in Lemma 2.4. 
Then the matrix function (3, which is defined in (2.1), satisfies the equality 



v{x) = \^'{x)j-i{x)*. 



(4.9) 




(4.10) 
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Proof. First, we fix an arbitrary I and rewrite (2.21) (for x < I) m the form 

j{x) = {E[^, Im,]){x). (4.11) 

It follows that 

i{EAE-'E<i>[){x) = ji{x) - 72(x)$i(+0). (4.12) 

Recalling that 7i(0) = and taking into account (2.17) and Remark 2.3, we 
have $i(+0) = 0. Therefore we rewrite (4.12) as 

i{EAE-^E^[) (x) = 7i(x). (4.13) 

Next, we substitute K = EAE'^ from (2.6) into (4.13) and ( using (2.4) and 
(2.8)) we get 

71 (x) = -^{x)j j\{tr(^E^[){t)dt. (4.14) 

Formulas (4.11) and (4.14) imply 

71 (x) = -7(x)j £ {E 7^,] ){tr {t)dt. (4.15) 

Because of (2.33), (2.34) and (4.15) we see that 

^{x)je{xr ^ 0, e{x) : = [V 0] + {E^'^{tr{E [$i /^,] ){t)dt 

(4.16) 

= [V o]+2'(5-Vi)(t)*[$i(i) 

where Sx ■= PxSPx- We shall show that 9 = (5. 

In view of (4.11) and the second equality in (4.16), we have 

e'{x) = [E^'^{xy^{x). (4.17) 

Therefore, (2.3) leads us to 

/3{x)je'{xy = 0. (4.18) 

17 



Furthermore, compare (2.3) with the first equahty in (4.16) to see that 

e{x) = x{x)l3{x), (4.19) 

where x{x) is an mi x mi matrix function, which is boundedly differentiable 
on [0, l\. Now, equalities (4.18) and (4.19) and the first relations in (2.3) and 

(4.6) yield that V = (i.e., x is a constant). It follows from (2.1), (4.16) 
and (4.19) that x(0) = 7^^, and so x = 7^^, that is, 9 = /3. Thus, (4.10) is 
immediate from (4.16). □ 

Remark 4.4 Because of (4.5), (4.17) and equality /3 we see that the 
potential v can be recovered via the formula 

v{x) ^ (iE^[yxy. (4.20) 

Formally applied, formulas (2.34), (2.39) and (4.20) yield 

vix)^{iS;'^[)ixr, (4.21) 

though one needs a proper "pointwise" definition of matrix functions S~^<^'^ 
for (4.21) to hold. 

The last statement in this section is a Borg-Marchenko-type uniqueness the- 
orem, which follows from Theorems 3.7 and 4.1. 

Theorem 4.5 Let ip and (p he the Weyl functions of two Dirac systems on 
[0, oo) with the locally bounded potentials, which are denoted by v and v, 
respectively. Suppose that on some ray = c'^z (c G M, > 0) the 
equality 

\\ip{z) - ip{z)\\ = ©(e^'''^) {Qz oo) for all 0<r<l (4.22) 
holds. Then we have 

v{x)^v{x), 0<x<l. (4.23) 
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Proof. Since Weyl functions are non-expansive, we get 

\\e-'^"'' {ip{z) - ip{z)) II < cie^"'"', '^z>C2>0 (4.24) 

for some ci and C2, and the matrix function e"^^'^^ {^ip{z) — (p{z)) is bounded on 

the hne Qz = C2- Furthermore, formula (4.22) imphes that e'"^^"^^ {ip[z)~(p[z)^ 
is bounded on the ray '^z — d^z. Therefore, applying the Phragmen-Lindelof 
theorem in the angles generated by the line '^sz — C2 and the ray '^z — d^z 
{^z > C2), we see that 

lle-2'''^ {ip(z) - if(z)) II < C3, Qz>C2>0. (4.25) 

Functions associated with ip will be written with a hat (e.g., $1). Because of 
formula (3.9), its analog for (p, $1 and the inequality (4.25), we have 

e2'(^-'')^($i(x) - ^i{x))dx\\ < C4, ^z>C2> 0. (4.26) 

Clearly, the left-hand side of (4.26) is bounded in the semi-plane Qz < C2 
and tends to zero on some rays. Thus, we derive 

pr 

/ e2^(^-'^)^($i(a;) -$i(x))cix = 0, i.e., $i(x) = $i(x) (0 < x < r). 
Jo 

(4.27) 

Since (4.27) holds for all r < I, we obtain $i(a;) = ^i{x) for < x < I. In 
view of Theorem 4.1, the last identity implies (4.23). □ 

5 Weyl function and positivity of S 

In this section we discuss some sufficient conditions for a non-expansive ma- 
trix function (p, which is analytic in C+, to be a Weyl function of Dirac 
system on the semi- axis. For the case that v is a, scalar and Weyl functions 
(f are so called Nevanlinna functions (i.e., ^(f > 0), sufficient condition for 
</? to be a Weyl function can be given in terms of spectral function [28,29], 
which is connected with (p via Herglotz representation. Furthermore, a posi- 
tive operator S is also recovered from the spectral function (see [43, Chapters 
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4,10] and [35]). The invertibihty of the convolution operators, which is re- 
quired in [4,5,25], provides their positivity too, and the spectral problem is 
treated in this way. Here, we again formulate conditions on the m2 x mi non- 
expansive matrix functions in terms of S. To derive our sufficient conditions, 
we apply the procedure to recover Dirac system from its Weyl function (see 
Section 4). Recall that $i in Section 4 is the Fourier transform of (p, that is, 
it is given by (4.1). 

First, consider a useful procedure to recover 7 from (3 as mentioned in 
Remark 4.2. 

Proposition 5.1 Let a given m-yXm matrix function I3{x) (0 < x < Z) he 
houndedly differentiate and satisfy relations 

Pio) = [V 0] , P'jP* = 0. (5.1) 

Then there is a unique m2 x m matrix function 7, which is houndedly differ- 
entiahle and satisfies relations 

7(0) = [0 Im,] , ih* = 0, ^jP* = 0. (5.2) 
Moreover, this 7 is given by the formula 

7 = 727, 7 = [71 Im,] , 71 = ^2(^1*)-', (5.3) 

where det f3i{x) 7^ 0, det (/„i2 — 71(2^)71(2;)*) 7^ and 72 can he recovered via 
the linear differential system and initial condition helow: 

72 = 727^71* (^m2 - 7i7i)"\ 72(0) = Im,- (5.4) 

Proof. Because of (5.1) we have l3j/3* = Imi (and so det^i 7^ 0). On the 
other hand, formula (5.3) implies 7j73* = 0. Therefore, we get 

7j;5* = 0, 7i7* < 0. (5.5) 

In particular, we see that det {l,n2 — 7i7i) 7^ 0. According to (5.4) we 
have det 72 7^ 0. Furthermore, any 7 satisfying (5.2) admits representation 
7 — 727 with some houndedly differentiable 72 such that 72(0) = and 
det 72(2;) 7^ 0. Thus, it remains to rewrite 7'j7* in the equivalent form 

(727 +[7271 0])j7* = 0, (5.6) 
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which, in its turn, is equivalent to the first equahty in (5.4). Clearly (5.4) 
uniquely defines 72. □ 



Now, we formulate the main statement in this section. 

Theorem 5.2 Let an 1712 x mi matrix function ^p{z) be analytic and non- 
expansive in C+. Furthermore, let matrix function $i(x) and operators Si, 
which are given by (4.1) and (3.8), respectively, be such that $1 is boundedly 
differentiable on each finite interval [0, /] and satisfies equality ^i{x) = for 
x < 0, whereas operators Si are boundedly invertible for all < I < 00. 

Then (f is the Weyl function of some Dirac system on [0, 00). The op- 
erators Sf^ admit unique factorizations 

Sf^ = El^,E^,i, E^^i ^1+ [ E^{x,t) ■ dt e B{L^'{0,1)), (5.7) 

Jo 

where Eq>{x,t) is continuous with respect to x, t and does not depend on I, 
and the potential of the Dirac system is constructed via formula 

v{x) = (iE^^i^'^Yx)*, 0<x<l. (5.8) 

Note that formula (5.8) in Theorem 5.2 is similar to (4.20). To prove this 
theorem we need an auxiliary Proposition 2.1 from [15] on operator S: 

Proposition 5.3 [15] Let ^i{x) be an m^. x mi matrix function, which is 
boundedly differentiable on the interval [0, I]. Then the operator S, which is 
given by (3.8), satisfies the operator identity (2.18), where H := [$i $2] is 
expressed via formulas 

$,e5(C-^L^^(0, /)), ($i/)(a;) = $i(a;)/, ^^f^lmj- (5.9) 

Factorization results from [22, pp. 185-186] yield the following lemma (see 
[15, Theorem 4.2]). 

Lemma 5.4 Let matrix function ^i{x) and operators Si, which are expressed 
via $1 in (3.8), be such that^i is boundedly differentiable on each finite inter- 
val [0, /] and satisfies equality $i(0) = 0, whereas operators Si are boundedly 
invertible for all < I < 00. Then the operators Sf^ admit factorizations 
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(5.7), where E^{x,t) is continuous with respect to x, t and does not depend 
on I. Furthermore, all the factorizations (5.7) with continuous E^{x,t) are 
unique. 

Now, we consider our procedure to solve inverse problem. 

Lemma 5.5 Let the conditions of Lemma 5.4 hold. Then the matrix func- 
tions 

(3^{x) := [I^, 0] + £ {ty [$i(0 Im,] dt, (5.10) 

7$(x) := [$i(x) Im,]+ [ E^{x,t) [^i{t) Im,]dt (5.11) 

Jo 

are boundedly differentiable and satisfy conditions 

/g$(0) := /3$(+0) = [I^, 0] , I3'^jl3l = 0; (5.12) 
7<i.(0) = [0 /^,] , i'^nl = 0- 7*J/3S = 0. (5.13) 

Proof. Step 1. The first equalities in (5.12) and (5.13) arc immediate 
from (5.10) and (5.11), respectively. Furthermore, (5.11) is equivalent to the 
equalities 

-f^{x) = (^E^^i [$i /„,J ) (x), 0<x<l (for all Z < oo). (5.14) 

Next, fix any < Z < oo and recall that according to Proposition 5.3 the 
operator identity 

ASi - SiA* = iUjn* (5.15) 

holds, where IT is expressed via (5.9). Hence, taking into account (5.7) and 
(5.14), and turning around the proof of (5.15) in Lemma 2.4, we get 

E^AE^' - {E^'yA*El = i7$(x)i / 7$(t)* • dt, i.e., 

Jo 

E^AE^' = ij^{x)j [ 7,.(t)* ■ dt (for E^ = E^^i). (5.16) 
Jo 
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Introducing the resolvent kernel r$ by E^^ — I + r$(x, t) ■ dt, we rewrite 
(5.16) in the form of an equality for kernels: 

px px 

Im2+ / {E^{x,r) + r^{r,t))dr+ / / E^{x,r)drT^{^,t)d^ 
Jt Jt Ji 

= -7*(^)J7*(*)*- (5.17) 

In particular, formula (5.17) for the case that x — t implies 

l^{x)h^{x)* = -Im2- (5.18) 

In a way quite similar to the first part of the proof of Proposition 4.3, we use 
equalities (5.14), (5.16) and $i(0) = to derive 

l^{x)jM^y = (5.19) 

(compare with (4.16)). Because of (5.7), (5.10) and (5.14) we have 

P'^ix) = (x)*7$(x). (5.20) 

In view of (5.19) and (5.20), is boundedly differentiable and the last 
relation in (5.12) holds. 

Step 2. It remains to show that 7$ is boundedly differentiable and the 
identity 7$i7|, = holds. For that purpose note that /3$ satisfies conditions 
of Proposition 5.1, and so there is a boundedly differentiable matrix function 
7 such that 

7(0) = [0 im,] , i'jT = 0; 7J/3; = 0. (5.21) 

Formulas (5.12) and (5.21) yield 

/3^j/3l = Irm, W* = -Im„ (5.22) 

respectively. That is, the rows of ^$ (the rows of 7) are linearly independent. 
Hence, the last relations in (5.21) and (5.22) and formulas (5.18) and (5.19) 
imply that there is a unitary matrix function cu such that 

7$(a;) = u!{x)'-y{x), a;(a;)* = u}{x)~^. (5.23) 
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Moreover, formulas (5.12), (5.21) and (5.22) lead us to the relations 



uju J ^ij 



V 

V* 



u{x) :— 



.7(a;) 



Uju*j = Im-, m(0) = Im- 



(5.24) 
(5.25) 



According to (5.24) and (5.25) the matrix function u{x) is the normalized 
by (1.3) solution of Dirac system, where v is the bounded potential and the 
spectral parameter z equals zero. In other words, /3$ and 7 correspond to v 
via equalities (2.1), and we can apply the results of Section 2. Therefore, by 
Proposition 2.2 and Remark 2.8 there is an operator E of the form (2.39), 
such that 



'm2 ■ 



EA = ft{x)j [ 7(i)* • dtE, 72 = EI„ 
Jo 

It follows from (5.23) and (5.26) that 

ujEA = i'y^{x)j / 7$(t)* • dtuE, 'y2,<s> = cuEIm^, 



(5.26) 



(5.27) 



where u denotes also the operator of multiplication by the matrix function 
u. On the other hand, formulas (5.14) and (5.16) lead us to 



E^A — ^7$ 



Jo 



dt E^, 72,$ — E^Ij 



1712 ■ 



(5.28) 



It is easy to see that c.l.s.lJ^Qlm (A* 1^2) = ^"^^[0,1). Hence, equalities 
(5.27) and (5.28) imply Eq> = uE. Taking into account that by (2.39) and 
(5.7) expressions for both operators E and E^ include a term /, we see that 



uj{x) = Im2, E^^E. 



(5.29) 



Furthermore, because of (5.23) and (5.29) we have, 7$ = 7, and so 7$ is 
boundedly differentiable and satisfies (5.13). □ 

Proof of Theorem 5.2. By the assumptions of Theorem 5.2 the conditions 
of Lemmas 5.4, 5.5 are fulfilled. Therefore, the theorem's statements about 
operators Si and E^ i are true. Furthermore, because of (5.18), (5.20), (5.24) 
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and equality 7$ = 7 we have v — v for v given by (5.8) and v from Lemma 
5.5. It follows that v is bounded on [0, /]. 7$ = 7 = 7 and operator E, which 
is recovered from v in Section 2, satisfies (5.29). Thus, we derived: 

V = v, 7$ = 7, E^ = E. (5.30) 

Using (2.17), (5.14) and (5.30) we see that $1, which is recovered from v and 
E in Section 2, coincides with $1 in the statement of theorem. In view of 
Proposition 3.2 and Corollary 3.8 there is a unique Weyl function Lpw of the 
Dirac system with the potential v and this Weyl function is given by (3.19). 

It remains to show that ipw equals the function (p, which generates $1 
via (4.1). RecaUing that (4.1) holds also for cpw: we see that 

/a 
-——{ip{^ + iri)-^w(^ + iri))d^ = 0, ri>0, (5.31) 

where l.i.m. stands for the entrywise limit in the norm of L'^{—r, r) ( < r < 
cxd). Therefore, we get = (f, that is, is the Weyl function of the Dirac 
system, where the potential is given by (5.8). □ 
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